The bivalve Thyasira cf. gouldi hosts chemoautotrophic symbiont populations with strain level diversity

Invertebrates from various marine habitats form nutritional symbioses with chemosynthetic bacteria. In chemosynthetic symbioses, both the mode of symbiont transmission and the site of bacterial housing can affect the composition of the symbiont population. Vertically transmitted symbionts, as well as those hosted intracellularly, are more likely to form clonal populations within their host. Conversely, symbiont populations that are environmentally acquired and extracellular may be more likely to be heterogeneous/mixed within host individuals, as observed in some mytilid bivalves. The symbionts of thyasirid bivalves are also extracellular, but limited 16S rRNA sequencing data suggest that thyasirid individuals contain uniform symbiont populations. In a recent study, Thyasira cf. gouldi individuals from Bonne Bay, Newfoundland, Canada were found to host one of three 16S rRNA phylotypes of sulfur-oxidizing gammaproteobacteria, suggesting environmental acquisition of symbionts and some degree of site-specificity. Here, we use Sanger sequencing of both 16S RNA and the more variable ribulose-1,5-bisphosphate carboxylase (RuBisCO) PCR products to further examine Thyasira cf. gouldi symbiont diversity at the scale of host individuals, as well as to elucidate any temporal or spatial patterns in symbiont diversity within Bonne Bay, and relationships with host OTU or size. We obtained symbiont 16S rRNA and RuBisCO Form II sequences from 54 and 50 host individuals, respectively, during nine sampling trips to three locations over four years. Analyses uncovered the same three closely related 16S rRNA phylotypes obtained previously, as well as three divergent RuBisCO phylotypes; these were found in various pair combinations within host individuals, suggesting incidents of horizontal gene transfer during symbiont evolution. While we found no temporal patterns in phylotype distribution or relationships with host OTU or size, some spatial effects were noted, with some phylotypes only found within particular sampling sites. The sequencing also revealed symbiont populations within individual hosts that appeared to be a mixture of different phylotypes, based on multiple base callings at divergent sites. This work provides further evidence that Thyasira cf. gouldi acquires its symbionts from the environment, and supports the theory that hosts can harbour symbiont populations consisting of multiple, closely related bacterial phylotypes.